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Abstract 



Solid material in a protoplanetary nebula is subject to vigorous redistribution 
processes relative to the nebula gas. Meter-sized particles drift rapidly inwards 
ne£ir the nebula midplane, and material evaporates when the particles cross a 
condensation/evaporation boundary. The material cannot be removed as fast in 
its vapor form as it is being supplied in solid form, so its concentration increases 
locally by a large factor (more than an order of magnitude under nominal con- 
ditions). As time goes on, the vapor phase enhancement propagates for long 
distances inside the evaporation boundary (potentially all the way in to the 
stcir). Meanwhile, material is enhanced in its solid form over a characteristic 
lengthscale outside the evaporation boundary. This effect is applicable to any 
condensible (water, silicates, etc.). Three distinct radial enhancement/depletion 
regimes can be discerned by use of a simple model. Meteoritics applications in- 
clude oxygen fugacity and isotopic veiriations, as well as isotopic homogenization 
in silicates. Planeteiry system applications include more robust enhancement of 
solids in Jupiter's core formation region than previously suggested. Astrophysi- 
cal applications include differential, time-dependent enhancement of vapor phase 
CO and H2O in the terrestrial planet regions of actively accreting protoplanetary 
disks. 
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1 Introduction 



Matter doesn't simply condense from a cooling protoplanetary nebula at its cosmic relative 
abundance and remain in place. Significant inward radial transport of solids occurs relative 
to nebula hydrogen (MorfiU and Volk 1984, Stepinski and Valageas 1997), and trace vapors 
migrate outwards to condensation fronts (Stevenson and Lunine 1988). Each compound has 
its own condensation front; water's is often called the snowline. Previous work has stressed 
the role that condensation of outwardly diffused vapor plays in enhancing the density of 
solids at the snowline (Stevenson and Lunine 1988). However, inward particle drift can also 
enhance the abundance of a vapor inside the condensation/evaporation boundary (Cuzzi et 
al 2003). Hence, given turbulent mixing, the density of solids just outside the boundary 
also grows. That particle drift can have this effect has been recognized (Morfill and Volk 
1984), but the magnitude of the effect has not been recognized. Here we show that, in 
general, inward particle drift is more effective than outward vapor diffusion at enhancing 
solids outside the condensation/evaporation boundary, and that in many cases of interest 
vapor is also strongly enhanced interior to the condensation/evaporation boundary. To 
emphasize the difference between this new process and previous work, we will often refer 
to the condensation/evaporation boundary as the evaporation front, which more accurately 
captures the directionality of the process described here. 

In this paper we present a minimal model that illustrates the key physical arguments. We 
divide a condensible solid into three size classes distinguished by their transport properties: 
vapor and small grains that are tightly coupled to the movements of the gas; large bodies 
that orbit unaffected by gas drag; and mid-sized particles (boulders, or rubble — typically on 
the order of a meter) that are strongly affected by both gravity and gas drag (Weidenschilling 
1977). The latter can drift orders of magnitude faster than the nebula gas and carry a net 
flux greatly exceeding that which is coupled to the gas. The particles evaporate when they 
reach the evaporation front, enhancing the abundance of the condensible in the vapor phase 
(Cuzzi ct al 2003). While the physics is applicable to volatilcs in general, we focus here on 
water as a volatile of special interest. The enhanced water vapor abundance spreads radially 
inwards from the evaporation front on a timescale which is short compared to the lifetime 
of the nebula, potentially determining the mineralogy of primitive meteorites. Planetesimal 
growth just outside the condensation/evaporation boundary, or snowline, provides a sink that 
ultimately depletes the vapor at all locations inside the boundary (Stevenson and Lunine 
1988). Ice enhancement outside the snowline can influence the location and timescales of 
giant planet core formation (MorflU and Volk 1984, Stevenson and Lunine 1988). Our model 
illustrates all these regimes of behavior. 

1.1 Nebula evolution and turbulence: 

Protoplanetary disks evolve; both their surface mass density ag and mass accretion rate 
M = 2nRagVn (where Vn is the nebula gas advection velocity) decrease with time over a 
period of several Myr (Calvet et al 2000). The disk is heated by gravitational energy release 
{GM/R) and illumination from the star. In recent models of actively accreting disks (Bell 
et al, 1997, Stepinski 1998; cf. also Woolum and Cassen 1999), midplane temperatures 
are hot enough in the terrestrial region to vaporize common silicates at early times, and 
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to vaporize water in the jovian region over a more extended period. The physical cause of 
nebula evolution remains problematic; turbulent viscosity is now thought to face difficulties 
(Stone et al 2000). However, turbulence can exist, providing diffusivity, without necessarily 
providing the viscosity needed to evolve the disk (Prinn 1990; Cuzzi et al 2001). Since 
diffusivity, rather than viscosity, is of prime interest here, we will assume nebulae which are 
weakly turbulent. As discussed below and by Cuzzi and Weidenschilling (2004), turbulence 
plays several roles: it diffuses grains and vapor down concentration gradients, often against 
the flow of nebular drift (Cuzzi et al 2003); it frustrates the growth of particles beyond 
about a meter in size; and it determines the midplane particle density (and thus particle and 
planetesimal growth rates). 

We presume a weakly turbulent nebula with effective viscosity vt — acH, where c is the 
sound speed, H is the scale height, and the parameter a is defined by the evolutionary mass 
accretion rate of the nebula: a = M /STrcTgcH . Hence the advection velocity Vn = 3acH/2R. 
Observations and models typically suggest that a ~ 10~^ — 10~^ and H ^ R/20. The 
turbulent diffusivity V = ut/ Prt is related to the viscosity by the turbulent Prandtl number, 
Prt- It is usually presumed that Prt = 1 in a turbulent nebula, but in general Prt ^ 1; e.g. 
Prinn (1990) suggests that Prt <^ 1 and that mixing is efficient even if the nebula evolves 
slowly. The characteristic velocity of large eddies is ~ c{a / PrtY^^ . Nebular evolution is 
driven by pt, while V describes mixing. 

1.2 Turbulence and particle growth to meter-size: 

Particle growth is easy up to meter-size, but problematic beyond. With or without turbu- 
lence, the relative velocities between sub-m-size particles are low and the first aggregates 
probably grow by simple sticking into porous, dissipative structures (Weidenschilling and 
Cuzzi 1993; Cuzzi and WeidenschiUing 2004). Most growth occurs as large particles sweep 
up smaller ones (Weidenschilling 1997). Under nonturbulent conditions, large particles sink 
into a high density midplane layer in which relative velocities are low, and subsequent growth 
to planetesimal sizes is rapid (Weidenschilling 1997, Cuzzi et al 1993). In turbulence, be- 
cause meter-sized particles couple to the largest eddies and achieve random velocities on the 
order of Vg, they remain in a layer of finite thickness /il ~ Vg/Qx ~ H(a/PrtY^'^. Even 
weak turbulence {a > 10~^, Prt = 1) keeps the density of the midplane layer and the ensu- 
ing particle growth rate substantially below their nonturbulent values (Dubrulle et al 1995). 
Furthermore, meter-sized particles collide with each other at speeds comparable to - me- 
ters to tens of m/sec - probably fragmenting into their smaller constituents. While particles 
approaching meter-size, having impact strength of 10^ erg/cm^ (Sirono and Greenberg 2000), 
can survive mutual collisions in turbulence of a = 10~^ (Cuzzi and Weidenschilling 2004), 
one suspects that further incremental growth may stall at the m-size limit, at least while 
turbulence this large persists. Other physics, however, may come into play {eg., Cuzzi et al 
2001). Detailed models of incremental growth, with realistic sticking and erosion based on 
laboratory experiments, tend to quickly produce broad power law size distributions which 
contain equal mass per decade of particle radius (Weidenschilling 1997, 2000). Assuming 
that growth beyond the meter-size range is frustrated as described above, and conserva- 
tively assuming ten decades of particle size (microns to meters would be six decades), we 
estimate that for extended periods of time, meter-sized particles have surface mass density 
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(Tl ~ O.lasol where asoi ~ 10 is the total surface mass density of sohds at some location. 
We will define (Jl/oSoi = /l ~ 0.1 as a key element of the model described in section 2. 

1.3 Radial drift: loss or transformation? 

The nebula gas, in general, has an outward pressure gradient, which counteracts solar gravity 
to a small degree; the ratio of these two forces is r] ^ H"^ / 2(? j^V]^ 2 x 10~^, where 
Vk is the local Keplerian velocity (Weidenschilling 1977, Nakagawa et al 1986, Cuzzi et 
al 1993) and 7 = 1.4 is the ratio of specific heats. The gas orbits more slowly than the 
solids at any location, and the ensuing headwind on particles causes them to drift radially 
inwards at velocities which depend on their size. Even weak turbulence ensures that local 
particle densities are too low to affect their drift velocities (Nakagawa et al 1986, Cuzzi and 
Weidenschilling 2004). Typical radial drift velocities are shown in figure 1; they are strongly 
dependent on particle size, but only weakly dependent on distance R from the Sun for nebula 
models such as that adopted [og = 1700(1 AU/R) g cm~^). Meter-sized particles experience 
the full headwind and are the most rapidly drifting; smaller particles experience a smaller 
headwind, and larger particles have increasing mass per unit area (Weidenschilling 1977). 
For comparison we show the range of gas advection velocity Vn at 5 AU and a = 10~^, which 
characterizes a range of radial disk density profiles. The drift velocity Vl ~ tjVk of m-size 
particles is orders of magnitude larger than Vn- It is usually inferred that such drifters are 
"lost into the sun" on fairly short timescales; however, this is not necessarily their fate — as 
we describe below. 

When particles drift across the location where the midplane temperature exceeds the 
sublimation temperature of one of their constituent species, that species will evaporate locally 
within distance AR ~ {m/m)VL where m is the particle mass and to its evaporation rate. 
For water, AR < 1 AU (Supulver and Lin 2000). Cyr et al (1998) have also modeled 
evaporation of drifting water particles, finding that they drift considerable distances before 
evaporating. However, these results are incompatible with those of Supulver and Lin (2000) 
and with our own estimates using the same vapor pressure expressions as Cyr et al. Even 
after discussions with J. Lunine (personal communication, 2003) we cannot determine the 
cause of this discrepancy. For silicates, AR is probably larger (Cuzzi et al 2003). This 
evaporation front eHect can produce significant enhancement of material. We have considered 
only volatiles which represent significant fractions of the total condensible mass at their 
evaporation front: common iron-magnesium silicates, at about 1400K (Cuzzi et al 2003) and 
water ice, at about 160K (this paper). The process is sketched in Figure 2 and described 
further in section 2. 

1.4 Caveats on assumptions 

The key parameter is uncertain. First, as discussed in section 1.2, and most recently 
by Cuzzi and WeidenschiUing (2004), we believe that the actual combination of nebula 
turbulence a and particle strength (Sirono and Greenberg 2000) allows particles to grow to 
meter-size. If this is not true, the appropriate value of Vl will decrease in proportion to 
the particle size-density product (Cuzzi and Hogan 2003). This would decrease the mass 
flux in "large" particles. Secondly, one of our simplifying assumptions was that growth was 
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log particle radius (cm) 

Figure 1: Radial drift velocity as a function of particle size, at a number of nebula locations 
differing by a factor of 1.5 (several labeled, in AU from the sun). The drift velocity at the 
peak is ~ tiVk- For comparison, a range is shown for the nebula advection velocity at 
it! =5AU with a — 10~^ (the range is associated with different nebula density profiles; cf. 
Cuzzi et al. 2003). 
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Evaporation Front 



Figure 2: Sketch illustrating inwardly drifting (blue) volatile material crossing its evaporation 
front Revi with midplane temperature Te^,. The surface density of large midplane solids is 
= fhCoO'g] the large inward drift flux of this material, aiYi, can't be offset by vapor 
removal processes CagVn + VagdC / dR until the concentration of the blue vapor C is much 
greater than nominal solar Co- The more refractory (red) material, here shown as a minor 
constituent, simply goes on drifting and growing. 
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truncated above meter-size, and that planetesimals don't grow. This led to our estimate 
of fi - the mass fraction in meter-sized particles - as simply the inverse of the number of 
decades of size in the "rubble" population between microns and meters. If turbulence is 
vanishingly small and planetesimals do grow, this simple logic loses its appeal. However, 
even in a planetesimal growth regime, coUisional erosion or breakup will continue to occur. 
In this regime, might be regarded as the coUisionally generated, small-size end of a mass 
distribution with equal mass per decade, extending to 1000km. 



2 Model for evaporation fronts 

Here we will construct a simplified 1-D model that illustrates the basic principles that govern 
transport of a condensible in the solar nebula. Let C represent the mobile fraction of a 
condensible species in the nebula. We will define this as the column density of the species 
(vapor or solid) divided by the column density of the nebula (chiefiy hydrogen and helium 
gas). We exclude from C material that has condensed on large planetesimals; large bodies 
will be treated as a stationary sink. 

We further subdivide C between the mass fraction in fast-drifting m-size rubble (sec- 
tion 1.2) and the complementary fraction 1 — in small grains (or vapor) that are strongly 
tied to the motions of the gas. In general we expect that both coUisional growth and dis- 
ruption of particles are faster processes than advection, so that the relative proportions of 
small grains to m-sizcd rubble is roughly constant where solids are stable. The evaporation 
front is defined hy R = R^y] inside the evaporation front /x, = 0. For reasons discussed in 
sections 1.2 and 1.4, we will assume that /l = 0.1 when R > R^v Transport of C is then 
described by 

where the inward radial mass fiux $ is the sum of nebula advection, diffusion, and particle 
drift, respectively: 

$ = (1 - h) Ca,V„ + ^^{(1 - h) Ca,} + hCa.VL. (2) 

In equation 1 the term on the right hand side represents the accretion sink onto large plan- 
etesimals; Tacc represents an accretion time. For reasons discussed above (see also Cuzzi and 
WeidenschiUing 2004), the m-size particles are much more concentrated toward the midplane 
than smaller size ranges which contain equal amounts of mass. We have therefore assumed 
that accretion onto large bodies is dominated by m-size particles. 



2.1 Steady state solutions and a likely transient case 

We simplify equations 1 and 2 by assuming steady state, constant coefficients, and cartesian 
geometry. The latter two assumptions introduce quantitative errors on the order of unity 
provided that nebular properties (other than those associated with the evaporation front) 
vary smoothly on the scale of i?, as is usually assumed in discussions of the nebula. We are 
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left with 



where as before = for i? < 

We then further simphfy equation 3 by placing all accretion onto planetesimals at the 
condensation front. This captures the spirit of the snowline without introducing a complete 
model of planetary accretion. With this simplification, equation 3 is directly solved analyti- 
cally. There arc four boundary conditions. At large distances R ^ Rev the nominal cosmic 
abundance is C = Co', i.e. C{R oo) = Co- At small distances R <C Rev there is no source 
of C. This precludes the purely mathematical solution in which outward diffusion balances 
inward advection for R < R^v Consequently C{R < Rev) = C, a constant. The other two 
BCs apply at Rev We assume that C is continuous across Rev, and we apply a flux jump 
condition across Rp 
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/f Cn rl n 
" = CVnCag, (4) 

Rev — 



where £ is a dimensionless sink factor integrated over the narrow band of planetesimals just 
outside of Rev, defined to make the sink term similar in form to other terms in equation (3). 
The steady state solution that results is 

C{R < Rev) — ECo 

C{R > Rev) ^Co[l + {E- 1) e{-'=(«-«-)}] (5) 

where 

, _ (1 - /l) K. + hVL 

{l-h)D 

and 

„ _ {l-fL)V„ + hVL 

^ = (i + £)v; ■ 

The factor E is the enhancement over cosmic abundance. If /l ^ 1, 

J, ^ (1 + hVL/Vn) Eo 

where the factor Eo is that of Cuzzi ct al. (2003). 

In steady state the entire region interior to Rev is enhanced over cosmic abundance (in 
the vapor) by the factor E. The distance scale 1/A; is closely related to E. It is like a skin 
depth. It represents the distance scale beyond Rev in which solids are enhanced. Note that 
if Pvt < 1, the skin depth deepens accordingly. 

Enhancements can be large. Using Vl ~ tjVk we can estimate that, for {f^, C) <C 1, 

E^Eo^lp^^^^^, (9) 
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which, for /l — 0.1 and 10~^ < a < 10~^, is a very large factor indeed. By contrast, Morfill 
and Volk (1984) got much smaller vapor phase enhancements (never exceeding solar) because 
they assumed particles which drift only at about the same rate their nebula was advecting 
{Vl = Vn), and because of their choice of outer boundary condition (their equation B7). 

We discern three regimes of interest for C {R, t) , shown schematically in figure 3. Regimes 
2 and 3 are the steady state solutions described by equation 5 above. In regime 2 (£ ^ 1), the 
entire region inwards of R^y is enhanced by E over solar. Regime 3 occurs when planetesimal 
growth is significant and a sink appears at R^^ {C > 1). If £ is big enough, the inner nebula 
can become depleted, essentially the result of Stevenson and Lunine (1988). Regime 1 
(sketched only conceptually in figure 3) is transient, because a certain amount of time is 
needed to reach steady state. At first, evaporated material is found only within a radial 
band of width AR (for water AR <C R, Supulver and Lin 2000). This transient solution 
propagates toward the star and approaches a steady state only after a time t^s ~ Rev/Vn ~ 
l/{3nari) fti iO/a orbit periods. Of course, tss will also depend on Pr^. For Rev=5 AU, 
Prt = 1, and a — 10~^ to 10~^, tgg ~ 0.5 - 5 Myr — long enough to be interesting for the 
chemistry of the early inner nebula (Cuzzi et al 2003). Depending on the rate at which £ 
grows, the nebiila might evolve from regime 1 through regime 2 into regime 3, or directly 
from regime 1 into regime 3. 

2.2 Global constraints on the model 

Naturally, the steady state enhancement regime can't persist for the entire duration of disk 
accretion. For example, in regime 2, with ~ 100, as much water is accreting onto the sun 
as hydrogen! This enhanced stage is limited in duration and intensity by (a) growth of the 
planetesimal sink at Rev, leading to emergence of regime 3, and (b) depletion of the ultimate 
source of the enhancement - outer solar system solids. 

In most nebula models, the surface mass density decreases as 1/R; thus if the nebula 
extended only to 50 AU, 10 times further than R^.^, and if all the solids in that region 
were to be carried into the region interior to Rev, only an enhancement factor of 10 could 
be achieved; with < 1, the hmit could be even lower. However, the true radial extent 
and mass distribution in the actual nebula are unknown; many protoplanetary disks are 
not tens, but hundreds of AU across. Furthermore, some nebula models {eg. Ruden and 
Pollack 1991) show the nebula surface mass density increasing outwards, due to the effects of 
radially varying viscosity. In such a case, the same global constraint allows enhancement by 
a factor of 400, even if the nebula only extended to 50 AU. The hkely time-variable nature of 
more realistic solutions should be kept in mind. Even while global source constraints limit 
the steady state solutions, large Eo might prevail over limited times and radial distances. 
Improving astronomical observations of the radial extent and surface mass distribution of 
protoplanetary nebulae will be helpful in establishing such global constraints. Overall, we 
do not feel that values of Eg 10 — 100 are unreasonable (especially in regime 1) but Eg 
could be smaller (especially in regime 2) because of global constraints. 
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Figure 3: Schematic of the radial (and temporal) variation of enhancement C /Co for "water" 
with an evaporation boundary at i?eD=5AU, taking for illustration i?o=20. In regimes 1 and 
2, there is no sink at Rev (>C = 0); regime 1 (dotted; schematic only) represents the transient 
situation, where the inner nebula retains C/Co = 1 for typically 40/q; orbit periods. Regime 
2 (dashed) is the steady state solution for £ = 0. As time proceeds and planetesimals grow 
in the enhanced solid density outside R^v, C increases; regime 3 (solid) illustrates the steady 
state solution for Eo=2Q and C = 100. 
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2.3 The sink: planet esimal growth outside the condensation / 
evaporation boundary 

The sink term L removes solid material from further radial evolution by accreting it onto the 
surfaces of immobile planetesimals just outside R^.^. Thus the local mass density of potential 
planet- forming objects increases. A detailed study of this process is well beyond the scope 
of this paper, but the following simple expressions illustrate the possibilities. 
The mass lost to the planetesimal sink can be written as 

MpL = 'I'KR^^ECoagVnC = MECoC. (10) 

Independently we can write the mass accreted by a narrow belt of Np^ planetesimals with 
radius tpl and total mass MpL as 

MpL = Npj^i^rl, {^-^^^^^ Al^e, (11) 

where ~ V^, ~ viVk is the relative velocity of sweep up, ^ is a sticking coefficient, and 
where we have ignored gravitational focusing (appropriate for tpl < 30 km). Random 

velocities for mctcr-sizcd particles are comparable to Vg, so hi/H ^ {ajPrtY^"^ (Cuzzi and 
Hogan 2003; Cuzzi and Weidenschilling 2004). Setting equations 10 and 11 equal to each 
other, some algebra leads to 



MpL \ UlPt 



1/2 



In evaluating equation 12 we assume /l = 0.1, ps — 1, Prt — 1, and Re,v — 5 AU. Smaller 
bodies, which present a larger surface area for a given mass, are more efficient sinks provided 
that they are large enough to be immobile (greater than 100 meters or so; the size is itself 
a-dependent). Without detailed accretion modeling, it is hard to go further. However, at 
this level of description, interesting ranges of values for £ and MpL can be estimated (Table 

1)- 

Equation 10 may be rewritten to estimate the planetesimal belt growth time (in orbit 
periods) : 

Mp^ Mp^ _ 2.4 iMpL\ ^^^(Mpj\nj^\ ^^g^ 



Mpi &T:^Rl^ECoagr}aC EaC \ , 

Equation 13 assumes a snowhne at R^v — 5 AU, where Og — 300 g/cm^ and Cg — 0.01. We 
used E ~ Eo/ (1 + C) for Jl ^ 1, and equation (9) with fi = 0.1. 

In the C <^ 1 regime, equation 13 may be combined with equation 12 to obtain a 
characteristic growth time for the belt (in orbit periods, shown in line 3 of table 1): 



i^iC « 1) ^ 2 X 10^ a'/'r' ( ^) (14) 
/for, ^ Vlkm/ ^ ^ 



MpL 



The transition from regime 2 to regime 3 can happen very quickly once nebula turbulence 
dies down {a decreases) near R^y] however, the poorly understood sticking coefficient ^ enters 
into all these timescales. 
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Table 1: Interesting ranges of planetesimal belt mass and loss factor 



a -5- 


io-» 


10-^ 


lo--* 


10-3 


10-2 


MpL/Mg, (for £ = 1) 


5 X 10-3 


0.15 


5 


150 


5 X 103 


£ (for MpL = Mffi) 


200 


7 


0.2 


7 X 10-3 


2 X 10-'' 


Mp^jMpL, periods 






200 


6 X 103 


2 X lO'^ 



* indicates the high-/", regime, wliere equation (14) is inappiicable 



The dependence of the planetesimal accretion sink on nebula a is illustrated 
in this table. In all cases, rpi=l km and ^ = 0.1 are assumed. The first and 
second rows are both solutions of equation 12. Row 1 shows that the mass of 
a planetesimal belt capable of providing £ = 1 ranges from very small values 
for low a to impossibly large values for large a (a minimum mass nebula at 5 
AU contains about 1 of solids in a band of radial width H) . Row 2 shows 
the values of £ provided by a belt with Mpl = . Row 3 is the solution to 
equation 14; the mass doubling time for a planetesimal belt in the low-£ regime 
is independent of the mass of the belt, and also shorter for lower o. For small 
a, accretion onto planetesimals is faster and the inner nebula quickly becomes 
more depleted. 
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3 Applications 



3.1 Meteoritics and astronomical observations: 

Nebulae which arc hot enough to evaporate sihcates near the midplane in the terrestrial 
planet region (M ~ lO^^M0/yr, or age ~ 10^ years; Bell et al 1997) arc probably also young 
enough that the drifting solids are more primitive and carbon-rich than chondrites. Evap- 
oration of silicates in the presence of 20-30% carbon by number may lead to the formation 
of abundant CO, with interesting mineralogical and isotopic implications (Cuzzi et al 2003). 
For the duration of regime 1, this enhanced silicate and CO vapor plume near i?ei; (silicates) 
need not be accompanied by a similarly enhanced component of water, because the water 
stripped out of the drifting solids at i?et) (water) remains at radii ^ i?ei, (silicates) until steady 
state is achieved some tgs = 40/a orbit periods later (regime 2). The potential duration of 
this dry, CO-and-silicate-rich inner solar system regime 1 seems to be comfortably longer 
than the apparent duration of the CAI formation era in the inner solar system (Cuzzi et al 
2003), which plausibly ends when the inner nebula cools to below the evaporation temper- 
ature of common silicates - and thus long before the chondrule era that apparently occurs 
1-3 Myr later {eg., Amehn et al 2002; Russell et al. 2004). 

Before planetesimal growth at i?e,, (water) creates a sink, regimes 1 and 2 can provide an 
enhanced abundance of H2O relative to hydrogen over a wide range of locations. This may 
help explain several aspects of chondrite chemistry indicative of elevated oxygen abundance, 
such as high levels of FeO in matrix olivines in both ordinary and carbonaceous chondrites 
(Nagahara 1984, Scott et al 1984, 1989; Wood 1988). It has traditionally been argued 
that enhancement of nebula gas in silicates of chondritic composition can provide the high 
oxygen fugacity required for high-FeO silicates to form in the nebula (Palme and Fegley 
1990). Recently, however, Fedkin and Grossman (2004) have shown that chondritic silicates 
are ineffective in this regard because they provide too much S (which competes for Fe) in 
addition to their O. Our mechanism enriches the nebula gas in H2O alone, and might provide 
the needed oxygen fugacity without the sulfur complications. In another application, Ciesla 
et al (2003) have suggested that fine grained sihcates can be aqueously altered in the nebula 
gas by shock waves, if the nebula gas is enhanced in H2O by something like a factor of 100. 
This level of enhancement is achievable, even if perhaps only regionally or for a limited time, 
under circumstances described here. Furthermore, the enhancement is probably temporally 
and spatially variable, depending on how the nebula evolves between the three regimes we 
have identified. In addition, A. Krot (2003, private communication) has noted that this 
enhancement can affect the Oxygen isotopic ratios in primitive meteorite minerals in a time- 
variable way, if outer solar system ice has different 0-isotopic composition from inner solar 
system silicates (eg., Lyons and Young 2004; Yurimoto et al 2004; Krot et al. 2004). 

Some recent astronomical observations seem to show abundant CO in the terrestrial 
planet regions of vigorously accreting protoplanetary nebulae (Najita et al 2003). The pres- 
ence of abundant CO might be associated with the evaporation front of primitive silicate- 
carbon material discussed above. In at least one case, the water content in the inner nebula 
seems to be low relative to CO (Carr et al 2003). This could be the signature of regime 1, 
or perhaps a very early stage regime 3. Future observations of this type, perhaps at higher 
spatial resolution, might help us determine evolutionary timescales and connect the current 



14 



properties of external protoplanetary nebulae with the record of the accretion process in our 
own. 



3.2 Planetary formation 

Formation of Jupiter has long been associated with the concept of a snowhne (Wuchterl et al 
2000). The nominal scenario for the formation of an icy jovian core in less than a few Myr, 
while the nebula gas is still present, requires that the surface mass density of solids in the 
formation region exceed that of a minimum mass nebula by almost an order of magnitude 
(Lissauer 1987). One well-known proposal for this enhancement is the cold finger effect, in 
which the entire water content of the inner nebula is diffusively transported to the snowline 
and frozen out there (Stevenson and Lunine 1988). With the assumptions of a vigorously 
turbulent inner nebula, a narrow condensation annulus, and no leakage back into the inner 
solar system (questioned by Sears 1993), the cold finger effect leads to an enhancement of 
solids outside the snowhne by a factor of 6-25 in about 10^ years. This can be expressed as 
a mass flux of roughly 7rRl^Coag/{Rl^/T>) ~ nCgagV ~ a few times 10~^ Earth masses per 
year for Pr^ ~ 1. The ratio of the mass flux to i?e^ due to solids drifting from o-uter regions 
(this paper), to that of solar abundance vapor diffusing from the mner solar system is 

Unless the turbulent Prandtl number is very small, i.e. unless turbulent transport is much 
larger than viscous transport, and given the validity of our particle size distribution argu- 
ments, inward particle drift and vapor retention would seem to be the dominant source for 
enhancement of solids near Re,^. 

Looking somewhat further beyond the boundaries of this paper, we suspect that evapo- 
ration fronts of low-temperature volatiles might also have important implications for their 
enhancement in the gaseous envelope of Jupiter, a problem highlighted by Owen et al (1999) 
and Atreya et al (2003). Another possibility, associated with i?eD (silicates), is the isotopic 
homogenization of a large amount of the silicate material which ultimately ends up in me- 
teorite parent bodies. The gross isotopic homogeneity of meteoritic silicates has been a 
persistent puzzle, because few nebula models evaporate silicates throughout the asteroid 
formation region. It has also been suggested to us (J. Chambers, personal communication, 
2003) that the process, operating at (silicates), might help explain the mass distribution 
in the terrestrial planet formation region. Some of these applications will be addressed in 
future papers. 

Summary: We show that nebula constituents will be enhanced in the vicinity of their 
condensation/evaporation boundaries Rev, due to rapid inward drift of solid material in the 
form of meter-sized boulders, and slow subsequent removal of the ensuing vapor. This evap- 
oration front effect modifies not only the surface mass density of solids available just outside 
Rev (useful for planet building) , but also the chemistry and mineralogy of material which re- 
sides well inside Rev (of potential importance to meteoritics and gas giant atmospheres) . The 
enhancements can be one or even two orders of magnitude, and probably vary on timescales 
of a million years or so - perhaps also exhibiting significant radial variation during that time. 
Some of these properties might be observable by astronomical observations. 
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Table 1 : Interesting ranges of planetesimal belt mass and loss factor 



a 


10-6 


10-5 


10-^ 


10-3 


10-2 


Mpl/Ms, (for £ = 1) 


5 X 10--^ 


0.15 


5 


150 


5 X 10^ 


£ (for MpL = Me) 


200 


7 


0.2 


7 X 10-3 


2 X 10-^ 


Mpl/Mpl, periods 


* 


* 


200 


6 X 103 


2 X 10^ 



* indicates the high-£ regime, where equation (14) is inapplicable 



The dependence of the planetesimal accretion sink on nebula a is illustrated 
in this table. In all cases, rpL=l km and ^ = 0.1 are assumed. The first and 
second rows are both solutions of equation 12. Row 1 shows that the mass of 
a planetesimal belt capable of providing £ = 1 ranges from very small values 
for low a to impossibly large values for large a (a minimum mass nebula at 5 
AU contains about 1 of solids in a band of radial width H) . Row 2 shows 
the values of £ provided by a belt with Mpl = Mq. Row 3 is the solution to 
equation 14; the mass doubling time for a planetesimal belt in the low-£ regime 
is independent of the mass of the belt, and also shorter for lower a. For small 
a, accretion onto planetesimals is faster and the inner nebula quickly becomes 
more depleted. 
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